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THE COST EFFECTIVENESS OF MINICOMPUTERS VS. 
MAIN FRAMES FOR STRUCTURAL ANALYSIS PROBLEMS 


Raphael T. Haftka 
Anneliese K. Von Mayrhauser 

Illinois Institute of Technology 
Chicago, Illinois 60616 


ABSTRACT 


A study of the cost effectiveness of mini¬ 
computers vs. main frames for structural analy¬ 
sis programs is described. The study compares 
the performance of several finite element pro¬ 
grams including SAP IV and SPAR. Most of the 
runs were performed with the Illinois Institute 
of Technology PRIME 400 minicomputer and the 
United Computing System UNIVAC 1100/81 main 
frame. Other computers were used selectively. 
The matrix of structural problems included 
beam, plate and shell problems and static, dy- 
pamic and non 
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MINICOMPUTERS PROMISE CHEAPER, more widely avail¬ 
able computing facilities, but they pose many 
problems, particularly to those with large calcu¬ 
lations in mind. The smaller main memory means 
that users must make more use of disk I/O. Many 
minicomputers have a smaller word size (8 or 16 
bits per word) with a devastating effect on accu¬ 
racy. (Storaasli and Foster (1)* report 4 digit 
accuracy on a PRIME 400 for a medium-sized prob¬ 
lem as compared to 8-13 digit accuracy on a COC 
CYBER 173.) Longer elapsed times leave more room 
for machine hardware failures during a run. At 
Berkeley, Pearson reports some theoretical calcu¬ 
lations required twenty-four hours or more, but 
the mean time between failures (MTBF) was also 
approximately twenty-four hours. Lack of ade¬ 
quate software libraries, debugging facilities, 
documentation and operations staff often hampers 
users of minicomputers. 

Several past studies of the performance of 
minicomputers and main frames have concentrated 
on the machine cost aspect. Chemists have taken 
the lead In-doing large scale calculations on 
small machines, perhaps because of their famili¬ 
arity with the use of minicomputers for labora¬ 
tory process control. Peter Lykos of the 1111- 


81 7 


nois Institute of Technology organized an Ameri¬ 
can Chemical Society Symposium on Minicomputers 
and Large Scale Computations in June, 1977 (2). 
Four of the papers from this meeting are partic¬ 
ularly relevant to the present study since they 
involve large minicomputers with a cost in 
$100,000 - $200,000 range: Pearson et. al. (3); 
Norbeck and Certain (4); Wagner et. al. (5); and 
Freuler and Petrie (6). These papers have shown 
that based on machine cost alone minicomputers 
may be two to four times more cost effective than 
main frames. However, the National Resource for 
Computation in Chemistry Committee has concluded 
that "Minicomputers alone cannot meet the needs 
of the chemistry community" ( (7), p. 1). 

The structural analysis community is also 
experiencing a shift from main frames to mini¬ 
computers (e.g., Swanson (8) ). For example, 
most of the major finite element packages are 
now available on the PRIME 400. 

The users of structural analysis software 
on minicomputers have conducted several studies 
that point to the economic advantages of using 
such machines (e.g., Storaasli and Foster (1), 
Conaway (9) ). These are based on data gathered 
by running structural analysis programs for vari¬ 
ous problems. The scope of the comparison was 
limited by the relatively small amount of data 
that was generated. 

The purpose of the present study is to 
extend the previous studies in several direc¬ 
tions. First, to use a wider mix of programs 
and problems for the comparison. Second, to 
include in the study several machines so that 
the results do not reflect only the charging 
algorithm at a particular Installation. Last, 
the present study also attempts to further assess 
the human factors involved in such a comparison. 

CHOICE OF STRUCTURAL ANALYSIS PROGRAMS AND PROB¬ 
LEMS 


Since the overwhelming majority of struc¬ 
tural analysts use finite element methods, the 
present study is limited to such programs. 

It is impossible to get a very accurate com¬ 
parison of the performance of the mini and main 
frame systems without an exhaustive battery of 
test problems. However, we believe that by judi¬ 
cious choice of the type of problems and the type 
of structural analyses employed, one can obtain a 
reasonably reliable comparison of performance. 

•Numbers in parentheses designate References at 
end of paper 
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The following are the elements in such a compari¬ 
son. 

TYPE OF ANALYSIS-The type of structural 
analysis to be included in the tests must include 
the ones which are most commonly in use. This 
study includes the following: 

(i) Linear static solution for displace¬ 
ment and stresses 

(ii) Linear eigenvalue analysis - calcu¬ 
lation of buckling loads and vibration modes 

(Hi) Nonlinear response due to large 
deformations and material nonlinearities 

TYPE OF PROBLEM-Three structural problems 
are used: 

(1) A simple cantilever beam 

(ii) A plate with a hole 

(1ii) A stiffened cylinder 

These problems are solved for the different 
analysis types. Several models are used for 
each problem ranging from a very crude model to 
a refined one. The number of degrees of freedom 
ranging from a few dozen for a crude model to 
more than a thousand for the refined one. 

COMPUTER PROGRAMS-Three computer programs 
are used; these programs are: 

(i) SAP IV - A general purpose finite 
element program developed at the University of 
California at Berkeley is probably the most 
widely used "free" (it costs $200.00) finite 
element code. 

(ii) SPAR - A general purpose finite ele¬ 
ment program developed by W. D. Whetstone which 
serves as a prototype of a commercially devel¬ 
oped code. 

(iii) TWODEL - A special purpose finite 
element program developed by 0. Malkus at I IT 
for large deformation analysis of 2 dimensional 
elasticity problems. It is a representative of 
in house codes. 

THE COMPUTER PROGRAMS 

SPAR - SPAR (10) Is a general purpose finite 
element program developed by W. 0. Whetstone for 
NASA, first at Lockheed and then at his own EISI 
company. The program has linear static analysis, 
eigensolutions for vibration and buckling and 
model response capabilities. It does not have 
nonlinear analysis or direct Integration capa¬ 
bilities. The program has a public version 
which is distributed by the government through 
COSMIC. It also has a proprietary version 
called EAL (for Engineering Analysis Language). 
The public versions for the UNIVAC and CDC sys¬ 
tems are maintained by the developer while the 
minicomputer versions (PRIME and VAX) are main¬ 
tained by NASA. The proprietary version is 
available on all four systems but only In its 
executable version. Because of the expense of 
acquiring the proprietary version and the diffi¬ 
culty of Instrumenting It without access to the 
source code, the public version was used for 
this study. 

SPAR Is a modular system composed of more 
than 20 small programs called processors. The 
processors communicate through a data base sys¬ 
tem which is also directly accessible to the 
user. (See Giles and Haftka (11) for more 


information.) 

The public version of SPAR was installed on 
the UNIVAC without any problem. Installation 
on the PRIME was much more troublesome. Diffi¬ 
culties were due to some bugs in the PRIME ver¬ 
sion that had to be corrected and due to the vir¬ 
tual memory system of the PRIME that did not seem 
to work well for very large arrays. Addition¬ 
ally, the PRIME company issues new releases of 
the operating system quite often. Many times 
the older version of SPAR did not work with the 
newer operating system and the programs had to 
be recompiled and reloaded (a non trivial effort 
because SPAR is composed of so many individual 
programs). 

SAP IV - SAP IV (12) is also a general pur¬ 
pose finite element program that has static, 
vibration and dynamic analysis capabilities. It 
has been developed by Wilson and his students at 
Berkeley and is available at nominal cost ($200.) 
to the public. It is probably the most widely 
used "free" finite element program. There are 
more advanced versions of SAP denoted as SAP V, 
SAP VI, etc. which are available at considerable 
cost ($9,000.00) from the University of Southern 
California. In the present study SAP IV is used. 

The program was originally written for a 
CDC system. However, it has been converted to 
other systems. On the UCS UNIVAC 1100/81, there 
are three versions of SAP IV. However only one 
of these is working. The program can be com¬ 
piled using UNIVAC's FORTRAN V compiler. An 
attempt to use the more efficient ASCII FORTRAN 
compiler was unsuccessful. A substantial change 
in the program and I/O format statements may be 
required to make SAP IV suitable for this com¬ 
piler. There also is an absolute version of 
SAP IV which is supported by UCS and is avail¬ 
able for problem solution. 

The PRIME version of SAP IV was generated 
by Feeser of RPI. The first version of SAP IV 
which we received from the PRIME users' library 
was highly mutilated. Another tape was then 
obtained from RPI, Troy, after several months' 
wait. This tape had about 60 lines of code 
missing towards the end of the STRETR SUBROUTINE. 
Two following routines were also destroyed. 
Luckily this particular piece of code was correct 
on the earlier tape, and we were able to patch 
the code to make it work. Both times we received 
800 BPI tapes. The PRIME installation at IIT 
has only a 1600 tape drive, however. The con¬ 
version was another non trivial task. 

TWODEL - TWODEL (13) is a special purpose 
finite element program for two dimehsional finite 
elasticity developed by Malkus at IIT. It is a 
relatively small program and was developed simul¬ 
taneously for the UNIVAC and PRIME systems. The 
problem description is built into the program so 
that refining the mesh necessitates changes to 
program patches. 

COMPUTERS AND CHARGING ALGORITHMS 

A study on the cost effectiveness of a mini¬ 
computer and a main frame can be based on a 
detailed economic study of the costs of owning 
and operating each machine. Such a study is 
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feasible for a particular type of machine and a 
particular owner. It is, however, difficult to 
generalize to the broad category of minicomputers 
and main frames. The approach taken in the 
present work is to use the charging algorithms in 
several installations as the measure of the cost 
of running any of the selected problems. In 
most cases the charging algorithms are devised to 
recover the total cost of owning and operating 
the system, and in some cases to show profit. 

Most of the runs were performed at the com¬ 
puter center at the Illinois Institute of Tech¬ 
nology. The school operates its own PRIME 400 
minicomputer and buys computing services in a 
UNIVAC 1100/81 from the United Computing System 
Corporation (UCS). A few runs were made also on 
the CDC CYBER 176 owned by a large manufacturing 
company (LMC) and on the UCS CDC CYBER 176. The 
performance of the programs on the first three 
computers was used to predict performance on a 
few additional systems. The following is a brief 
description of the computer systems that were 
used in this study, their charging algorithms and 
the method used to predict performance on them. 

IIT - PRIME 400 - The PRIME 400 is a 
medium size minicomputer which is used at IIT to 
support interactive computing. The present con¬ 
figuration has 1.25 Megabytes of memory and two 
disk drives. The system is very busy from about 
9 AM Monday through Saturday with about 20 to 30 
users. However, most of the users are not CPU 
intensive. As a result the CPU Intensive struc¬ 
tural analysis programs often show very good 
response times even during the day. In running 
on the PRIME 400 it was found that the response 
time and 1/0 times are quite sensitive to the 
number of users while the CPU time Is not. To 
account for the variability runs were performed 
for several user environments on the PRIME and 
the performance averaged. 

The charging algorithms for the system is 
given in Table 1 (considerable discounts are 
available at non-business hours). While most of 
the users do not pay the charges, some do. The 
charging algorithm is intended to cover the total 
cost of owning and operating the system and pro¬ 
viding user support. Total billings are about 
$16,000.00 a month. 

NASA - LANGLEY RESEARCH CENTER PRIME 400 - 
The NASA Langley Computer Center operates sev¬ 
eral minicomputers Including a PRIME 400 com¬ 
puter. The charging algorithm (see Table 1) 
was devised to recover the capital and oper¬ 
ating costs. The NASA PRIME 400 supports a 
small number of users; however, the users run 
more CPU intensive jobs than the I IT users. 

Based on these considerations it was 
assumed that program performance Is similar on 
the ITT and NASA systems. The cost of running a 
job on the NASA PRIME was calculated based on 
this assumption. 

NASA - LANGLEY RESEARCH CENTER CYBER 173 - 
The NASA Langley Research Center operates sev¬ 
eral CDC main frames including two CYBER 173 
computers which are used to support interactive 
computation. The charging algorithm Is devised 
to recover capital and operating expenses. 

One of us (Haftka) has an extensive experi¬ 


ence using the SPAR program on both the I IT and 
NASA systems. Based on this experience the 
following assumptions were used to predict cost 
on the NASA CYBER 173. 

1. The CPU time on the CYBER 173 is 1.5 
times that of the CPU time on the UNIVAC 1100/81. 

2. The 1/0 charges are calculated based on 
the number of reads and writes from the SPAR 
data base available in the SPAR output. It was 
assumed that these numbers are similar for both 
systems, and that the charge per disk access is 
1.1 times the minimum charge per access. This is 
based on the feature of the charging algorithm 
which is relatively insensitive to the number of 
words in each read or write operation. 

3. The amount of core required is assumed 
to vary from 70K octal for the smallest number 

of nodes to 12QK for the largest number of nodes. 

4. Based on the above considerations the 
cost of running SPAR on the NASA CYBER 173 is 

Cost = (1 + 0.048A)(0.0136T 1 + O.OO3ni0) 

where 

T 1 = UNIVAC 1108/81 CPU time (sec) 

niB = combined number of disk reads and 
writes reported by SPAR 

A « core storage (In units of 10K octal 
words) 

LARGE MANUFACTURING COMPANY CDC CYBER 176 - 
A few runs of SPAR were made on a CYBER 176 which 
is owned by a large manufacturing company. The 
charging algorithm is given in Table 1. It 
appears to be fairly high on CPU and low on 1/0 
and core storage charges. In predicting costs 
on this machine it was assumed that CPU times can 
be predicted from the UNIVAC 1100/81 results 
based on published data (14) rating the UNIVAC 
1100/81 at 1800 K0PS and the CYBER 176 at 
9300 KOPS. It was assumed that the ratio of I/O 
time to CPU time is the same on both machines. 
Because the I/O cost on the CYBER 176 is low, 
even a large error in this assumption is not 
expected to change the cost much. For SAP IV 
the data for the UCS CYBER were used for the 
LMC CYBER estimates. 

CONCORDIA COLLEGE COMPUTER NETWORK (CCCN) - 
CCCN provides administrative, instructional and 
research computing services to a large number of 
educational institutions. The service Includes 
a large software development staff and the 
charges are computed to recover costs of operat¬ 
ing the'center. Interactive and batch services 
are provided by a UNIVAC 90/80 main frame 
Installed in December, 1 980. The UNIVAC 90/80 
is roughly equivalent to an IBM 370/158 in raw 
processing power. The prediction of costs of 
running on the CCCN UNIVAC 90/80 were based on 
the charging algorithm given in Table 1 with the 
following assumptions. 

1. CPU times are 2.25 longer than on the 
UNIVAC 1100/81. This Is based on published data 
(14) rating the UNIVAC 1100/81 at 1800 KOPS and 
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the UNIVAC 90/80 at 900 KOPS. 

2. Core requirements were estimated at 
200 K bytes for the small problems and at 300 K 
bytes for the large problems. 

IIT RESEARCH INSTITUTE VAX 11/780 - The 
I IT Research Institute (IITRI) operates a OEC VAX 
11/780 with 3 megabytes of core memory and 650 
megabyte disk space. The IITRI computer is not 
intended for uses involving much "number 
crunching" so that it does not have a floating 
point processor. It is not heavily loaded so 
that time is available to outside users. The 
charging algorithm is given in Table 1. 

RESULTS 

BEAM PROBLEM - A cantilever beam was mod¬ 
eled by plane beam elements (E24 elements in 
SPAR) and three vibration modes and frequencies 
were calculated using the SPAR program. The num¬ 
ber of nodes was varied from 5 to 600 and with 
three degrees of freedom per node; the maximum 
number of degrees of freedom is 1800. However, 
even for 600 nodes the problem is not costly to 
run because of the very small band width associ¬ 
ated with a one dimensional problem. 

For small problem sizes (up to 25) the prob¬ 
lem can be run in single precision. For larger 
number of nodes the limited double precision op¬ 
tion in SPAR (double precision used only for 
assembling the stiffness matrix) must be employed. 

This simple problem exposed a bug in the 
PRIME version of SPAR. The program could not 
calculate more than two vibration modes and fre¬ 
quencies. The problem occurred in the subspace 
iteration method and seemed to indicate that the 
initial vectors generated by a random number rou¬ 
tine were not linearly independent. On the 
UNIVAC there was no problem to get the required 
three lowest frequencies. 

The CPU, I/O, response times and cost for 
the beam problem are given for the UNIVAC 1100/81 
in Table 2 and for the IIT PRIME 400 in Table 3. 
The cost of running the problems is calculated 
based on the prime time rate. The total cost 
proved to be very close to a linear function of 
the number of nodes for both computers. It was 
therefore possible to predict accurately the cost 
of 600 node run by extrapolating the cost of the 
5-280 node runs. It is therefore assumed that 
the 5-600 node results can be used safely to 
extrapolate cost up to 1200 nodes. The actual 
and predicted cost results for both computers are 
also shown in Figure 1. From Tables 2 and 3 and 
Figure 1 it is dear that the beam runs are about 
twice as expensive to run on the.UCS UNIVAC 
1100/81 as on the IIT PRIME 400. The response 
time or turnaround time Is favorable on the PRIME 
for the small problems but the UNIVAC is quite 
faster for the larger problems. 

The performance of other computers for the 
beam problem was actually measured or estimated 
and the results are summarized in Table 4. It 
can be seen from Table 4 that the costs of run¬ 
ning the problem on the NASA PRIME 400 are com¬ 
parable to those at IIT. However, comparison 
with the other main frames reveal costs which are 
significantly lower than those of the UNIVAC 


1100/81. Part of the difference may be due to 
the fact that the COC CYBER machines are not 
operated for profit and do not have marketing 
expenses. The low cost on the CYBER 176 reflects 
low 1/0 charges on that machine compared to the 
NASA CYBER 173. On the CYBER 176 about $1.40 of 
the $2.24 cost is CPU related and the rest is 
1/0 related cost. On the NASA CYBER 173 only 
$0.57 of the $8.60 cost is CPU related and the 
rest is 1/0 related. In assessing the total cost 
results in Table 4 the main difference is not 
between minicomputer and main frames. Rather it 
is between the service bureau computers (the two 
UNIVAC machines) and the user owned machines 
which have lower charges. 

PLATE PROBLEM - A rectangular plate with 
a hole was modeled with quadrilateral plane 
stress finite elements. Because of symmetry only 
one quarter of the plate is modeled, see Fig¬ 
ure 2. The plate is subjected to uniform loads 
and the stresses in it are calculated for the 
cases where the number of nodes varies between 
35 and 594. The problem was run with the SPAR 
and SAP programs. CPU and I/O times response 
times and cost on the UNIVAC 1100/81 are given. 
From Tables 5 and 6 we see that the performance 
of the two programs is similar. The turnaround 
time or. the UNIVAC is much better than that of 
the PRIME. The cost on both machines are com¬ 
parable. 

It is possible to predict the cost of the 
594 node run from the first four runs by a quad¬ 
ratic polynomial with less than 10% error. It is 
expected, therefore, that it is possible to use 
the data in Tables 5 and 6 to predict with such 
a polynomial the cost of running problems up to 
about 1200 nodes. The cost data and the predic¬ 
tive curve are given for SPAR in Figure 3 and 
for SAP IV in Figure 4. It is seen from Fig¬ 
ures 3 and 4 that the cost of the plate problems 
Increases much more rapidly on the PRIME than on 
the UNIVAC. VAX results show a similar behavior. 

The results from Tables 5 and 6 were used 
to predict performance on other computers and 
the results are given in Table 7. The results in 
Table 7 do not seem to indicate a difference 
between the minicomputers and the main frames 
but rather that the cost on user owned COC 
machines is substantially lower than on the other 
computers. 

The plate problem was now analyzed for large 
deformations using the TWODEL program. The plate 
was assumed to be made of rubber modeled by the 
Mooney-Rivlin fi law with c, = 0.3 x 10 a psi and 
c- * 0.6 x 10° psi. The plate is stretched to 
twice of its original length and the solution 
performed with 10 Incremental steps. CPU and 
I/O times, response time and cost for the UNIVAC 
1100/81 are given in Table 8 and for the PRIME 
400 in Table 9. The TWODEL program uses only 
Incore operations so that the program size is 
very sensitive to the problem size. The required 
core size on the UNIVAC is also shown in Table 8. 

As can be seen from Figure 5 which compares 
the cost of running TWODEL for various problem 
sizes, the problem is much cheaper to run on the 
PRIME 400 than on the UNIVAC 1100/81. Further¬ 
more, the 289 node problem was the largest that 
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could be run on the UNIVAC 1100/81 because of 
core limitations. The PRIME 400 with its vir¬ 
tual memory could handle larger problems. The 
excellent performance of the program on the 
PRIME 400 is attributed to two factors. 

First, the large core requirements on the 
UNIVAC 1100/81 are very expensive (fifty per¬ 
cent of the cost of the 289 run was core cost). 
Second, the program was written by the developer 
on both machines. As a result it does not suffer 
from the deterioration in performance that 
afflicts programs that have been converted to 
different machines by users who are not as know¬ 
ledgeable about the program as its developer. 

This applies to the SPAR and SAP programs. A 
comparison of the CPU times on the PRIME 400 and 
the UNIVAC 1100/81 for the plate problem shows 
the following ratios. The PRIME 400 CPU times 
are 7-15 times longer for SAP IV, 7-13 times 
longer for SPAR, but only 5-8 times longer for 
TW0DEL. The core sizes given in Table 8, were 
used to predict performance on other computers. 
For C0C machines it is expected that there would 
be no need for double precision: The core 
requirements were accordingly slashed by 1.8 
and the same ratio was used for the CPU times. 

The comparison with other computers is given 
in Table 10. It is seen that the cost on the 
two UNIVAC computers is very high and the cost 
on the LMC CDC computer is very low compared to 
the PRIME 400. 

STIFFENED CYLINDER PROBLEM - A stiffened 
cylindrical shell, see Figure 6, modeled by 
plate and beam elements, because of syimetry 
only one half of a 90° segment of the shell was 
modeled. The lowest four vibration frequencies 
and the buckling Toad of the cylinder were cal¬ 
culated using the SPAR program. Three models 
were used with 5x5, 10x10 and 15x15 grids of 
elements. This corresponds to 36, 121 and 256 
nodes, respectively. The largest model has 
about 1500 degrees of freedom. 

The cylinder problem revealed two bugs In 
the SPAR program. The stress analysis which is 
a preliminary to the buckling calculation was 
not correct when combined membrane-bending ele¬ 
ments (type E43 in SPAR) were employed. As a 
result it was necessary to use twice as many 
plate elements; pure membrane elements (type 
E41) and pure bending elements (type £42). 
Unfortunately, the mass matrix was not calcu¬ 
lated properly with this replacement. As a 
result we had to perform the calculations in 
two separate runs, one using double elements for 
buckling analysis and one using combined ele¬ 
ments for the vibration analysis. This did not 
Increase the cost of the run substantially but 
was an Inconvenience for the analyst. 

The result of CPU and I/O times, response 
time and cost are given in Table 11 for the 
UCS UNIVAC 1100/81 and in Table 12 for the I IT 
PRIME 400. Comparison with other computers is 
given in Table 13 for the largest problem of 
256 nodes. It can be seen that for this problem 
the performance of the PRIME 400 is dismal and 
that of the COC very good. 

The experience with the same problem with 
the SAP IV program was quite different. The 


PRIME 400 version ran with no difficulties 
while the UNIVAC 1100/81 and the UCS CDC 
CYBER 7600 encountered convergence difficulties. 
These difficulties were bypassed by change of 
parameters for the 36 node problem but difficul¬ 
ties remained for the larger problems. It is 
possible that these could be resolved too, but 
it required a degree of familiarity with the 
SAP IV program that the authors do not have. 

The cylinder problem was next analyzed for 
transient response by direct integration on the 
SAP IV program. The response to an oscillating 
normal load was calculated for forty time steps 
(each of 0.0125 sec). Since the cyclinder prob¬ 
lem is rather expensive, it was not run for all 
sizes on all computers. For the 256 node prob¬ 
lem the amount of scratch disk space required 
by the program was of the order of five megabyte 
on the PRIME 400 and was impossible to run in 
periods when the disk utilization was high. On 
the UNIVAC the program automatically switches to 
use magnetic tapes in such a situation. This 
results In very large I/O times and high charges. 
The results are given in Table 14. This problem 
illustrates best how the problem size influences 
resource consumption and charges. That is, for 
the smallest problem overhead charges dominate 
which are lowest on the minicomputers and make 
them cheaper. For the medium sized problem CPU 
consumption is dominant and the main frame is 
cheaper, whereas for the large problem the addi¬ 
tional core requirements were accommodated more 
cheaply by the computers with virtual memory. 

CONCLUDING REMARKS 

A study of the performance of structural analysis 
programs on main frame and minicomputers has not 
demonstrated a clear cut advantage of either 
type of computer. Difficulties were encountered 
on both types of computers for different prob¬ 
lems. For the minicomputers some of these dif¬ 
ficulties arose from portability problems when 
structural analysis packages that were developed 
on main frame computers were converted imper¬ 
fectly to minicomputers. Another factor is that 
software is easier to get for a main frame. As 
the problem sizes increased, the minicomputers 
showed reduced efficiency for those programs 
which were not developed for them; but on the 
other hand, some of the main frames were unable 
to run those problems adequately, too. The mini¬ 
computer clearly outperformed the main frame for 
the one program which has been developed on the 
minicomputer. 

When costs are compared there was also not 
much difference between the main frames and the 
minicomputers except in the case of the user 
owned CDC CYBER machines. This probably 
reflects a highly competitive marketplace where 
the decrease in the cost of raw computation has 
stimulated an Increase in programming and com¬ 
puter services offered. Main frames no longer 
hold a monopoly on structural analysis soft¬ 
ware (or most other software for that matter) 
and have to price their services adequately to 
compete with the minicomputers. 

If one is merely interested in the cost, 
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the user owned machines which just charge to 
cover their costs naturally fare better than 
machines owned by organizations which need to 
make a profit and may have marketing and other 
related expenses. This reduces the problem to 
one of microeconomics and price theory (15). 
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Table 1 

Charging Algorithms (dollars) 






m 


i 




1. IIT PRIME 400 


2. NASA Langley Research Center - 
PRIME 400 

3. UCS UNIVAC 1100/81 


4. Concordia College UNIVAC 90/80 

5. Large Manufacturing Company CDC CYBER 176 


6. NASA Langley Research Center - 
CDC CYBER 173 

7. UCS CDC CYBER 176 


0.02 (T, + T 2 ) + T 3 
T.j * CPU time (sec) 

T 2 * I/O time (sec) 

T 3 * wall time (hours) 

30T-| + 15T 2 

T^ * CPU time (hours) 

T 2 = wall time (hours) 

0.18T 1 + 0.0011A(T 3 + T 2 ) 

T 1 * CPU time (sec) 

T 2 = I/O time (sec) 

A = (core/512) words 

(0.156 + 0.0000248A)T 
T * CPU time (sec) 

A = core (kilobytes) 

0.2T 1 + 0.0003A (T 1 + T ? ) + .03 * T 2 
T^ = CPU time (sec) 

T 2 = I/O time (sec) 

A = core (kilowords) 

See section on "Computers and Charging 
Algorithms," pages 2-4. 

Confidential 


8. IIT Research Institute VAX 11/780 


0.05T.J + 5T Z 

Tj = CPU time (sec) 

T 2 = connect time (hours) 


Table 2 


Beam Results on the UCS UNIVAC 1100/81 Using SPAR 


Number of 

CPU 

I/O Time 

Turnaround Time 

Cost 

Nodes 

(sec) 

(sec) 

(min) 

(dollars) 

5 

5.06 

15.77 

3:21 

1.83 

25 

5.74 

16.12 

2:50 

2.07 

60 

7.07 

18.17 

6:04 

2.66 

90 

8.30 

19.79 

2:41 

3.16 

125 

10.29 

22.81 

4:26 

4.07 

280 

16.70 

32.59 

7:17 

6.90 

450 

22.52 

41.73 

3:16 

9.53 

600 

30.73 

59.68 

7:45 

13.80 

1200 * 

60.76 

72.10 

8:55 

30.89 


* predicted 
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Table 3 


Beam Results on the IIT PRIME 400 (Averages) 


Number of 

CPU Time 

I/O Time 

Response 

Cost 

Nodes 

(sec) 

(sec) 

Time 

(dollars) 

5 

14.70 

17.68 

1:03 

0.67 

25 

19.33 

17.43 

1:03 

0.73 

60 

30.65 

28.45 

2:37 

1.2'’ 

90 

37.07 

22.84 

1:57 

1.23 

125 

47.82 

37.25 

6:02 

1.80 

280 

94.24 

63.11 

5:41 

3.25 

450 

141.07 

132.87 

11:04 

5.16 

600 

197.38 

125.76 

13:47 

6.68 

1200 * 

422.92 

262.02 

26:36 

14.14 


* predicted 


Table 4 


Comparison of the 600 and 1200 Node Beam Problems on Various Computers 
(1200 Node Results Are Given in Parentheses) 


Computer 

CPU 

I/O 


(sec) 

(sec) 

11T 

, 197.38 

125.76 

PRIME 400 

(422.92) 

(262.02) 

NASA LRC * 

197.38 

125.76 

PRIME 400 

(422.92) 

(262.02) 

UCS UNIVAC 

, 30.73 

59.68 

1100/81 

(60.76) 

(72.10) 

NASA LRC * 

46.1 


CYBER 173 

(91.14) 


LMC 

* 6.05 

24.5 

CYBER 176 

(11.76) 

(35.4) 

Concordia 

69.14 

-- 

College * 

UNIVAC 90/80 

(136.71) 



Turnaround Time Total 

(min) Cost (dollars) 


13:47 

6.68 

(26:36) 

(14.14) 

13:47 

5.08 

(26:36) 

(10.17) 

7:45 

13.80 

(8:55) 

(30.89) 


8.60 


(20.78) 

.. 

2.24 


(3.98) 


11.13 


(22.34) 


* estimated 
+ average 


Haftka and 
Von Mayrhauser 


8 








Plate Linear Analysis with SPAR on Various Computers 


Number of 

CPU 

1/0 

Turnaround 

Cost 

Nodes 

(sec) 

(sec) 

(min:sec) 

(dollars) 

UNIVAC 1100/81 

35 

4.80 

14.25 

1:27 

1.69 

72 

6.98 

15.55 

1:44 

2.46 

143 

11.89 

18.36 

2:08 

4.22 

285 

24.46 

25.27 

2:48 

8.69 

594 

60.58 

44.78 

4:36 

21.94 

1200 

170.27 

102.31 

9:23 

62.89 

IIT PRIME 400 

35 

32.11 

19.18 

1:23 

1.05 

72 

58.40 

29.91 

2:58 

1.81 

143 

123.14 

44.96 

4:30 

4.18 

285 

285.14 

123.65 

8:56 

8.39 

594 

805.06 

679.25 

41:24 

28.71 

1200 

2500.40 

1370.10 

144:15 

126.21 


Table 6 - Plate Linear Analysis with SAP IV on Various Computers 


Number of 

CPU 

1/0 

Turnaround 

Cost 

Nodes 

(sec) 

(sec) 

(min:sec) 

(dollars 

UNIVAC 1100/81 

35 

4.52 

1.11 

1:13 

0.90 

72 

6.77 

1.93 

1:51 

1.56 

143 

11.55 

3.54 

2:12 

2.97 

285 

25.38 

8.19 

3:08 

7.15 

594 

78.59 

26.94 

10:08 

23.61 

1200* 

316.58 

122.14 

18:54 

98.67 

IIT PRIME 400 

35 

30.82 

25.65 

7:55 

1.26 

72 

66.40 

41.15 

11:00 

2.33 

143 

166.12 

47.15 

21:27 

4.62 

285 

.396.80 

82.61 

25:32 

10.01 

594 

1186.52 

216.10 

54:08 

28.95 

1200* 

3791.00 

685.74 

112:12 

108.60 

UCS CYBER 176 

35 

.54 

1.88 

— 

0.96 

72 

1.00 

3.08 

-- 

1.92 

143 

1.69 

4.48 

-- 

3.36 

285 

4.10 

10.82 

-- 

7.20 

594 

11.96 

38.21 

-- 

21.60 

1200* 

32.96 

99.18 

-- 

70.81 

IIT Research Institute VAX 11/780 

33 nor 

+ ,♦* 

0:25 

0.52 

72 

20.94 

— 

0:30 

1.05 

143 

59.62 

— 

3:30 

2.98 

285 

182.43 

— 

5:28 

9.12 

594 

568.33 

— 

12:41 

28.42 


* estimates 
+ averages 

** jobs were run in batch mode, no connnect time charges 
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Table 7 


Comparison of the Cost of Linear Analysis of the 594 (1200) Node 


Plate Model 

on Various Computers (in Dollars) 



Computer 

SAP 

IV 


SPAR 

System 

594 

1200 

594 

1200 

I IT PRIME 400 + * 

28.95 

108.60 

28.71 

126.21 

NASA LRC PRIME 400 

23.38 

59.64 

17.06 

56.90 

UCS UNIVAC 1100/81* 

23.61 

98.67 

21.94 

62.89 

NASA LRC CYBER 173 

★★ 


6.84 

18.90 

LMC CYBER 176* 

4.23 

11.39 

2.85 

7.82 

Concordia College UNIVAC 90/80 

28.98 

116.42 

22.28 

69.57 

UCS CYBER 176 

21.60 

70.81 


-- 

IIT Research Institute VAX 11/780 

28.42 

★★ 

— 

— 


+ average of several runs 

* predicted results 

** data insufficient for prediction 


Table 8 


Plate Nonlinear Analysis on the UCS UNIVAC 1100/81 


Number of 

CPU Time 

1/0 Time 

Core Storage 

Turnaround 

Cost 

Nodes 

(sec) 

(sec) 

(Kwords) 

Time (min) 

(dollars) 

35 

19.97 

0.35 

22.9 

32:47 

4.61 

72 

54.15 

0.66 

28.4 

38:01 

13.13 

143 

156.91 

1.29 

42.5 

55:10 

42.68 

289 

549.69 

1.78 

76.0 

137:21 

189.33 

600* 

2030.59 

3.59 

175.7 

650:43 

1151.80 

* estimates 



Table 9 




Plate Nonlinear 

Analysis on the 

IIT PRIME 400 


Number of 

CPU 


I/O Time 

Response 

Cost 

Nodes 

Time 


(sec) 

Time (min) 

(dollars) 

35 

156.01 


2.74 

7:08 

3.29 

72 

379.95 


3.81 

10:09 

7.84 

143 

978.05 


19.47 

29:59 

20.45 

289 

2693.59 


69.70 

70:50 

56.44 

600* 

8693.00 


287.91 

204:14 

193.06 

* estimates 



Table 10 




Comparison of the Cost (in Dollars) of Nonlinear Analysis of the 
289 and 600 Node Plate Models on Various Computers Using the TWODEL Program 


Computer 

System 

CPU 

Time (sec) 

289 600* 

Core* 

289 

600* 

289 

Cost 

(dollars) 

600’ 

IIT PRIME 400 

2694. 

8693. 

_ 

_ 

56.44 

193.06 

NASA LRC PRIME 400* 

2694. 

8693. 

— 

-- 

40.13 

123.50 

UCS UNIVAC 1100/81 

550. 

2030.6 

76.0 

175.7 

189.33 

1177.66 

LMC CDC CYBER 176* 

58.8 

217.1 

42.2 

97.6 

12.52 

49.80 

Concordia College 

IlNTUflr* 

1238. 

4570.7 

320. 

740. 

203.03 

794.20 


UNIVAC* 


* estimates; + UCS UNIVAC 1100/81 and LMC CDC CYBER 176 core requirements are given in K words. Con¬ 
cordia College UNIVAC core requirements are given in K bytes. 
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Table 11 


Cylinder Buckling and Vibration Analysis on the UCS UNIVAC 1100/81 


Number of 

CPU Time 

1/0 Time 

Turnaround 

Cost 

Nodes 

(sec) 


Time (min) 

(dollars) 

36 

27.28 

19.34 

7:09 

7.36 

121 

110.65 

28.16 

10:11 

29.21 

256 

302.59 

48.45 

19:52 

79.97 



Table 12 




Cylinder Buckling 

and Vibration Analysis on 

the IIT PRIME 400 


Number of 

CPU Time 

1/0 Time 

Turnaround 

Cost 

Nodes 

(sec) 

(sec) 

Time (min) 

(dollars) 

36 

579.93 

343.28 

64 

19.53 

121 

3197.10 

1184.54 

156 

90.23 

256 

10784.0 

6589.98 

774 

360.38 



Table 13 




Comparison of the 256 Cylinder Problem on Various Computers 


Computer 

CPU Time 

I/O Time 

Total Cost 

l/l 

S 

1 — 

(sec) 

(dollars) 

I IT PRIME 400+ 

10784 

6590 

360.38 

NASA LRC PRIME 400* 

10784 

6590 

283.37 

UCS UNIVAC 1100/81 

302.59 

48.45 

79.97 

Concordia College 

680.8 

— 

110.96 

UNIVAC 90/80* 



12.64 

EMC CDC CYBER 176* 

58.2 

24.2 

NASA LRC CYBER 173* 

453.9 

-- 

29.55 




Table 14 




Direct Integration Transient Response of Stiffened Cylinder 


Number of 

CPU Time 

I/O Time 

Turnaround 

Cost 

Nodes 

(sec) 

(sec) 

Time (min) 

(dollars) 

36 - PRIME 400 

157.1 

19.1 

4 

3.59 

UNIVAC 1100/81 

12.3 

10.7 

4 

6.12 

CYBER 176 

7.8 

26.6 


14.40 

VAX 11/780 

65.12 

— 

3.1 

3.26 

121 - PRIME 400 

2795 

528.8 

147 

68.93 

UNIVAC 1100/81 

133.7 

66.1 

14.9 

42.50 

VAX 11/780 

1049.38 

— 

28.7 

52.47 

256 - PRIME 400 

10411 

1658 

470 

249.20 

UNIVAC 1100/81 

493.8 

2423.2 

120 

556.70 
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Fig. 2 - A typical model for plate with hole problem 
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Fig. 6 - A typical cylinder model 
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